Polaritonic characteristics of insulator and superfluid phases in a coupled-cavity array 



OO 

o 
o 



Oh' 
I 



> 

o 
o 



E. K. IrishQ C. D. Ogden, and M. S. Kim 
School of Mathematics and Physics, Queen's University Belfast, Belfast BT7 INN, United Kingdom 

(Dated: March 5, 2008) 

Recent studies of quantum phase transitions in coupled atom-cavity arrays have focused on the 
similarities between such systems and the Bose-Hubbard model. However, the bipartite nature of 
the atom-cavity systems that make up the array introduces some differences. In order to examine 
the unique features of the coupled-cavity system, the behavior of a simple two-site model is studied 
over a wide range of parameters. Four regions are identified, in which the ground state of the 
system may be classified as a polaritonic insulator, a photonic superfluid, an atomic insulator, or a 
polaritonic superfluid. 

PACS numbers: 73.43.Nq, 42.50.Pq, 03.67.-a 



I. INTRODUCTION 

Recently a family of models for coupled arrays of 
atom-cavity systems has attracted considerable atten- 
tion [E H, H, 0, [E]. Building on the success of opti- 
cal lattice experiments in demonstrating the quantum 
phase transition between superfluid and Mott insulator 
states [B], these proposals have been inspired by experi- 
mental advances in photonic crystals 0, optical micro- 
cavities and superconducting devices [1]. Theoreti- 
cally, they offer a fascinating combination of condensed 
matter physics and quantum optics. Much of the work 
so far has focused on the possibility of creating quantum 
phase transitions in systems that permit manipulation 
and measurement of individual lattice sites. On the prac- 
tical side, several applications in quantum information 
processing have already been proposed. These include 
generation of entanglement [lO|, cluster state quantum 
computation [111, Il2| | and transfer of a qubit through an 
array |l3|. 

The simplest version of the coupled-cavity model con- 
sists of a series of electromagnetic cavities, each contain- 
ing a single two- level system (qubit or atom), coupled in 
such a way that photons may hop between adjacent cav- 
ities [1, H, 01 ■ Another model under consideration uses 
four-level atoms in a configuration commonly exploited 
for electromagnetically induced transparency, with each 
cavity containing multiple atoms [1]. The case of several 
two- level atoms per cavity has been explored as well Islj. 
Studies have been carried out in the microscopic regime, 
with only a few cavities P, [1, 0] > as well as the thermo- 
dynamic limit, in which the number of cavities goes to 
infinity 

Most of the previous work on phase transitions in 
coupled-cavity systems has emphasized similarities to the 
Bose-Hubbard model Indeed, the four-level-atom 

model can be mapped exactly onto the Bose-Hubbard 
Hamiltonian a two-component Bose-Hubbard model 
has also been derived from the four-level-atom sys- 
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tem [15|. In other coupled-cavity models, evidence has 
been found for a quantum phase transition between Mott 
insulator and superfluid states, analogous to that in the 
Bose-Hubbard model i, i, i, [3 ■ 

By contrast, our goal in this paper is to identify some 
of the unique features of the coupled-cavity system. To 
that end we have chosen to take a microscopic approach, 
building up from the well-understood Jaynes-Cummings 
model. Specifically, we consider a system of two cavi- 
ties containing a total of two excitations. Previous work 
has demonstrated that many-body effects appear in finite 
systems of only a few cavities, including signatures of the 
superfluid-insulator phase transition 1, 3, JJ. This holds 
true even for a two-cavity system. (For this reason we 
will use some of the language of quantum phase transi- 
tions, particularly the terms "insulator" and "superfluid" 
for localized and delocalized states, respectively. How- 
ever, it should be understood that in our usage these 
terms refer to states of a small finite system, not true 
phases in the thermodynamic sense.) The advantage of 
this approach is that the dimension is small enough that 
exact numerical solutions are easily found and, perhaps 
more importantly, some analytical approximations can 
be used. In this way we hope to identify characteristics 
of the coupled-cavity system that can be further explored 
in larger systems as well as in the thermodynamic limit. 

One of the principal ways in which the coupled-cavity 
model differs from the Bose-Hubbard model is that two 
types of particles are involved. The effective on-site re- 
pulsion is provided not by a fixed classical potential but 
by the interaction of photons with the atom(s) in each 
site. This interaction depends not only on the strength 
of the atom-photon coupling but also on the detuning 
between their frequencies. The detuning, then, provides 
an additional parameter for the system. To further com- 
plicate the picture, the fundamental excitations in the 
absence of hopping are not bosons but rather entangled 
states of atoms and photons, known as polaritons. The 
relative weights of the atomic and photonic components 
of the polaritonic states change with the detuning. As a 
result, both the effective repulsion and the nature of the 
particles in the system depend on the detuning parame- 
ter. 
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In this paper we map out the parameter space of the 
two-cavity system over a wide range of values for the 
atom-cavity detuning and the photon hopping rate. The 
negative detuning and large hopping regimes, in partic- 
ular, have not been explored in depth previously. We 
identify an interesting insulator-superfluid transition in 
the limit of large negative detuning, which occurs at a 
hopping strength equal to the magnitude of the detun- 
ing. The character of this transition is distinctly different 
from the transition at small hopping and positive detun- 
ing studied by Angelakis et al. Q. 

Previous studies of coupled-cavity systems have uti- 
lized the variance of the total excitation number to dif- 
ferentiate between insulatorlike and superfluidlike states. 
We also examine the atomic excitation number variance, 
which allows us to distinguish between purely atomic or 
photonic ground states and states that are composed of 
polaritons. Using these two measures we identify four 
distinct types of states: polaritonic insulator, photonic 
superfluid, atomic insulator, and polaritonic superfluid. 
In the small hopping limit, the change from atomic insu- 
lator to polaritonic insulator to photonic superfluid oc- 
curs smoothly as the detuning is increased. The situation 
is quite different around the transition in the regime of 
large negative detuning: the transition from atomic in- 
sulator to polaritonic superfluid to photonic superfluid is 
discontinuous to lowest order in the atom-cavity coupling 
parameter. This difference highlights the complexity that 
can arise from the combination of individual atom-cavity 
dynamics and the coupling between the sites. 



II. COUPLED-CAVITY MODEL: SMALL 
HOPPING PICTURE 

In order to keep the complexity of the system to a 
manageable level, we consider the simplest possible case, 
consisting of just two identical cavities. Each cavity sup- 
ports a single field mode and contains a single two-level 
atom. Photons are allowed to hop between the two cavi- 
ties. The Hamiltonian for the two-cavity system is given 
by {h = 1) 

H = [wca]aj + UJa \e.j) {ej\ 

+ dia] \gj) {ej\ + a-j \ej) {gj\)] + A{a\a2 + 4ai), 

(1) 

where ujc and uja are the cavity and atom frequencies, re- 
spectively, g is the atom-cavity coupling strength, and A 
is the hopping strength. The operator dj (aj) is the low- 
ering (raising) operator for the field in cavity j. The 
states \gj) and jcj) represent the ground and excited 
states, respectively, of the atom in cavity j. Hence the 
operator \gj) {ej\ (|ej) {gj\) is the atomic lowering (rais- 
ing) operator for cavity j. 

The first term of Eq. ([T|) gives the internal Hamil- 
tonian for the atom-cavity systems. Individual cavi- 



ties are described by the Jaynes-Cummings model, using 
the rotating- wave approximation [TtI [isj . The first two 
terms in the sum correspond to the internal energies of 
the field and atom. The interaction between the atom 
and the field, given by the third term, contains only so- 
called "energy-conserving" terms, in which an excitation 
of the field (atom) is accompanied by a deexcitation of 
the atom (field). Finally, the last term of Eq. ([T]) de- 
scribes the hopping of photons between the two cavities. 

In the absence of hopping {A = 0) the eigenstates of 
the individual cavities are given by the polaritonic states 



|0.) - \g^) |0, 



(2) 



l^i ) = siuy |e,) \ in- l)j) -cosy \n,) , (3) 

9 
\nt) = coSy |e,) |(n- 1),) +siny Ig,) \n,) , (4) 

where i — 1,2 denotes the cavity number, \n) (n — 
1,2,3,...) is a photon number state, and tan0„ = 
2gy/n/ A, where A = coa — the detuning. The ener- 
gies of these states are given by 



Ang^ 



(5) 
(6) 



The total number of excitations in the system Af = 
d\di + 0,2(12 -I- |ei) (ei| -I- |e2) (62], is conserved. In this 
paper the analysis is restricted to the case of exactly two 
excitations. There are eight possible states for the sys- 
tem, divided into five subspaces where states in the same 
subspace have degenerate energies. In order of increas- 
ing energy, the subspaces are ® |l2^)}, {|2j~) ® IO2), 
|0i) ® 12-)}, ® |1+), |1+) ® II2-)}, {|2+) ^ IO2), 

|0i) (g) |2J)}, {|i:;^) (g) |lj)}. It is important to note that 
the ordering of the subspaces with respect to energy is 
always the same regardless of the parameter values. How- 
ever, the energy differences between the subspaces change 
significantly. Figure [1] illustrates the energy levels in the 
three limiting cases of zero detuning, large positive de- 
tuning, and large negative detuning. 



III. TOTAL EXCITATION NUMBER 
VARIANCE 

Having established the model, we next turn to the 
problem of selecting a measure (analogous to an order 
parameter) that can distinguish between superfluid-like 
and insulator-like states. For the Bose-Hubbard model in 
the mean-field limit the expectation value of the boson 
destruction operator is typically used as the order param- 
eter [19,]. In the Mott insulator state each site contains 
a fixed number of particles and the expectation value of 
the destruction operator vanishes, whereas in the super- 
fluid state the particle number per site is not fixed and 
thus the expectation value becomes nonzero. However, 
our system is restricted to exactly two excitations. In 
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FIG. 1: Energy levels for the two-cavity system in the absence of hopping {A = 0): (a) zero detuning (A — 0); (b) large positive 
detuning {A/g ^ 1); (c) large negative detuning {—A/g S> 1). 



this case the expectation value of any destruction oper- 
ator is identicaUy zero. Therefore we will first look at 
the "order parameter" utilized in Ref. 3, which is the 
variance of the total excitation number in a single cavity. 
This quantity neatly captures the essence of the transi- 
tion. In the insulator state the number of excitations per 
cavity is sharply defined and has zero variance. How- 
ever, in the superfluid state each cavity has a finite prob- 
ability of containing any number of excitations, resulting 
in a nonzero variance. The two cavities are, of course, 
completely equivalent; for definiteness we will work with 
cavity 1. 

The excitation number in cavity 1 and its variance are 
defined as, respectively, 



Ni = a\ai + |ei) (ei| 



(7) 

(8) 



A contour plot of AA^i as a function of the detuning A 
and the hopping A for the ground state of the system is 
shown in Fig. [21 

To begin with, we consider the transition studied by 
Angelakis et al. [3] , which occurs in the region defined by 
A > and A <^ g. At A = 0, the interaction between 
the atom and the field mode in a given cavity shifts the 
frequency of the cavity mode. This creates a photon 
blockade effect, prohibiting additional photons from en- 
tering the cavity [2l| . The photon blockade leads to 
a large energy gap between the lowest two subspaces in 
the two-cavity system [see Fig. [TJa)]. When the hopping 
is weak {A/g <^ 1), the ground state of the system is ap- 
proximately (X) as illustrated in Fig. [31 Contain- 
ing exactly one excitation per cavity, this state is anal- 
ogous to the Mott insulator state in the Bose-Hubbard 
model [Ti |. 

At A = the state is a fully entangled state of the 
atom and photon [Eq. ([3])]. The inset of Fig. [3] demon- 
strates that, in the atom-cavity basis, this state has equal 
probabilities for atomic and photonic excitations. There- 
fore the ground state of the system at zero detuning and 
small hopping may be described as a polaritonic Mott 
insulator state. 

As the detuning is increased, the energy gap becomes 
smaller and the photon blockade is destroyed. In the limit 




FIG. 2; Plot of AA^i in the ground state of the two-cavity, 
two- excitation system. Black corresponds to ATVi — (Mott 
insulator state) while white corresponds to AA^'i = 0.5 (su- 
perfluid state). Throughout the paper we have taken g/uja = 



A/g ^ 1 the two lowest-energy subspaces of the two- 
cavity system become degenerate in energy [Fig. [Hb)]. 
The lowest-order effect of the hopping is to lift the de- 
generacy, resulting in a unique ground state. This ground 
state consists of a superposition of polaritonic states such 
as that illustrated in Fig. [H However, the nature of the 
individual cavity eigenstates is also altered by the change 
in the detuning, as seen in the inset of Fig.[31 In the limit 
A/g 00 we have \n~) « — \g) \n) and the ground state 
becomes [51) ® [52) jli) <^ jb) - |(|2i) d) (Oa) + \0i) ® 
[22))]. This is a delocalized photon state, i.e. a photonic 
superfluid. A state of this form is obtained by apply- 
ing two iterations of the delocalized creation operator 
-y={a\ — £12) to the vacuum state [l^. This is exactly the 



V2^ 

ground state of the Bose-Hubbard model with two sites 
and two excitations, in the limit of large hopping. 
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FIG. 3: (Color online) Probability distribution, broken down 
by subspace, of the ground state of the system as determined 
by numerical diagonalization of the Hamiltonian. The inset 
shows the probability distribution among states with purely 
photonic, purely atomic, and mixed excitations. Parameter 
values are A = 0, A — g/100. This is a Mott insulator state 
of polaritons. 
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FIG. 5: (Color online) As Fig. [1 but with A = -Wg, A = 
g/lOO. This is an insulator state composed almost entirely of 
atomic excitations. 



site, reminiscent of a Mott insulator. 
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FIG. 4: (Color online) As Fig. [3 but with A = lOg, 
A = g/100. This is a superfluid state that is almost entirely 
photonic in nature. 



A completely different situation occurs when the de- 
tuning takes on large negative values. In this case the 
energy gap between the two lowest subspaces becomes 
larger rather than smaller, as seen in Fig.[ljc). The plot 
of ANi (Fig. [2]) indicates that the ground state of the 
system remains in a Mott insulator state. This state is 
shown in Fig. O Again, though, the nature of the atom- 
cavity states changes with the detuning. In the limit of 
large negative detuning, w |e) |0) (see the inset of 
Fig. [5]) . Thus the ground state changes from a polaritonic 
insulator state at A = to an atomic insulator state at 
— A/g ^ 1. The explanation for this is quite simple. 
From the definition of the detuning, A = Ua — iOc, it is 
evident that when the detuning is negative the energy of 
the atoms is smaller than that of the photons. Therefore, 
the state of minimum energy is that in which only the 
atoms are excited. Because each atom is restricted to a 
single excitation, the purely atomic state is a localized 
state containing a definite number of excitations on each 



IV. ATOMIC EXCITATION NUMBER 
VARIANCE 

In the analysis so far we have identified atomic and po- 
laritonic insulator regions as well as a photonic superfluid 
state. A natural question, then, is whether there exists a 
region in which the superfluid state exhibits polaritonic 
characteristics. In order to answer this question, we must 
first define a measure that quantifies the degree to which 
a state is polaritonic in nature. 

Within an isolated cavity, the polariton states are char- 
acterized by a combination of atomic and photonic exci- 
tations. The total number of excitations is conserved and 
thus has zero variance. The variances in the numbers of 
atomic and photonic excitations depend on the detun- 
ing, reaching a maximum at A = when the atomic and 
photonic degrees of freedom are maximally entangled and 
dropping off to zero in the large detuning limits where the 
states become either atomic or photonic in nature. 

In the single-cavity case, the atomic excitation number 
variance and the photon number variance behave simi- 
larly. However, in the coupled-cavity system, the photon 
number variance is nonzero in the photonic superfluid 
state as well as in polaritonlike states. For this reason 
the photon number variance is not particularly helpful 
for our purposes and we shall not consider it here. The 
atomic excitation number variance, on the other hand, is 
zero in both the atomic insulator state and the photonic 
superfluid state, and thus provides a useful measure of 
the polaritonic nature of the state of the system. 

A plot of the atomic excitation number variance ANai , 
where Nai = |ei) (ei|, is shown in Fig. [51 There are two 
regions in which ANai = 0. The first, with A < and 
A < —A, corresponds to the atomic insulator state in 
which both atoms are excited. The second, which has 
yl > — A, corresponds to the photonic superfluid state 
in which both atoms are in the ground state. When 
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FIG. 6: Plot of AA'^Ai in the ground state of the two-cavity, 
two-excitation system. Black corresponds to AA^ai = 0, while 
white corresponds to AA''i — 0.25. 



A <^ the atomic excitation number variance of the 
ground state reduces to that of the polariton in the 
Jaynes-Cummings model. The atomic excitation number 
variance peaks around A = and drops off rapidly as the 
magnitude of the detuning is increased. As A increases, 
the height of the peak in AA'^i remains roughly con- 
stant, but the position of the peak follows the boundary 
of the insulator-superfluid transition A w —A. 

The atomic excitation number variance identifies re- 
gions of polariton-like behavior, but it does not distin- 
guish between insulator and superfluid states. In order 
to isolate the polaritonlike superfluid region, we take the 
product of the total excitation number variance AiVi, 
which is nonzero in the superfluid state, and the atomic 
excitation number variance AiV^i, which is nonzero in 
states with polaritonic characteristics. The result is 
shown in Fig. [71 

It is apparent from Fig. [7] that the superfluid region 
does indeed overlap to some extent with the region of po- 
laritonlike behavior. This region may be identified as a 
superfluid state that is, to some degree, polaritonic in na- 
ture. One such state is shown in Fig[51 All five polariton 
subspaces are occupied, indicating superfiuid behavior. 
The inset demonstrates that photonic and atomic exci- 
tations coexist, consistent with the idea that the particles 
involved in the superfluid state are polaritonic. 
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FIG. 7: Plot of the product of the total and atomic exci- 
tation number variances AA'^iAA''^! in the ground state of 
the two-cavity, two- excitation system. Black corresponds to 
ANiANai = 0, while white corresponds to AA^iAA^ai « 0.1. 
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FIG. 8: (Color online) As Fig. [3] but with A = -lOg, A = 
lOg. This represents a superfluid state with strong polaritonic 
characteristics. 



V. COUPLED-CAVITY MODEL: SMALL 
INTERACTION PICTURE 

The approach taken in Sec.|lTl in which A is treated as 
a small parameter and the system is described in terms of 
polariton states, accounts for most of the behavior of the 
coupled-cavity system. The atomic and polaritonic in- 
sulator states and the photonic superfiuid state all arise 
at small hopping {A <C g). The shift between atomic 
and polaritonic behavior of the insulator state is com- 
pletely determined by the atom-field interaction within 
each cavity, described by the Jaynes-Cummings model. 
In the limit of large positive detuning, the lowest-energy 
Jaynes-Cummings eigenstates when ^ = are photonic 
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in nature, and the superfluid state arises from the zeroth- 
order perturbative effect of the hopping. 

The polariton-hke superfluid state, on the other hand, 
appears only when A > g. In this case the description 
of Sec. |TT] breaks down. A different approach, in which 
g rather than A is taken as the small parameter, yields 
greater insight into the large hopping regime and the ap- 
pearance of the polaritonic superfluid. 

The Hamiltonian ^ may be split into three parts: a 
cavity Hamiltonian He, consisting of the harmonic oscil- 
lator term for each cavity plus the photon hopping term; 
the atomic Hamiltonian Ha', and the atom-cavity inter- 
action Hamiltonian Hi. These are given by 

He — ^ Wcajflj -I- A(a\a2 + alai), (9) 

3 = 1-2 

^ Wa|e,)(e,|, (10) 

Hi= 9{a] \gj) {ej\ + % \ej) {gj]). (--^-^^ 

Similarly, the basis states may be divided into three 
groups. The states that contain only photonic excita- 
tions are 

IV'cl) = |5ili)®l52l2), (12) 

|V'c2) = |5l2l) ® |ff202) , (13) 

iV'ca) = I5i0i) ® I5222) . (14) 

Only one state contains atomic excitations alone: 

\i^a) = \eiOi) ® \e2O2) . (15) 

Finally, there are four states that each contain one pho- 
tonic excitation and one atomic excitation, which are 

= leili)® |.g202), (16) 

1^.2) = I51O1)® |e2l2), (17) 

\i^^3) = leiOi) ® |.g2l2) , (18) 

|^^4) = |5ili)'» |e202). (19) 

We consider first the case g — and look for the ground 
state of the system. Hc + Ha is block diagonal in the basis 
given by the states p^ - p^ and may be diagonalized 
exactly. The resulting eigenenergies are {2ujc — 2A, 2ujc, 
2lUc + '2A, 2wc-H2A, 2lUc~A+A, 2lUc--A+A, 2uJc + A+A, 
2ujc + A + A}. 

Noting that A > and —00 < A < +00, three differ- 
ent regimes may be identified. These regimes are distin- 
guished by the relative values of A and A. 

When A < —A, the ground state energy is 2ujc ~ 2A, 
corresponding to the atomic insulator state |V'a)- It may 
be seen, therefore, that the atomic insulator state is not 
restricted to ^ <C 5. As discussed earlier, the negative de- 
tuning provides an energy gap between the ground state 
and first excited state. The hopping A must be on the 
order of |A| in order to overcome the gap and create a 
superfluid state. 
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FIG. 9: (Color online) Probability distribution in the atom- 
cavity basis of the numerically determined ground state (left- 
hand bars) compared with that given by Eq. (|2H) (right-hand 
bars). Parameter values are A = —lOg and A = lOg. 

On the other hand, when A > —A, the ground state 
energy is 2wc — 2A. This corresponds to the eigenstate 
IV'ci) = 7f IV'ci) - 5(|V'c2) + iV'cs)), which is the photonic 
superfluid state. Evidently, the photonic superfluid state 
is not found only in the large positive detuning regime: 
it emerges when the hopping becomes large enough, re- 
gardless of the value of A. 

The situation becomes slightly more complicated when 
A = —A. There are four degenerate eigenstates, two of 
which are IV'ci) |V'a)- The other two are given by 
W.i) = ^(IV'^2) - IV-a)) and \^',,) = ^(l^a) - IV'.a)). 
In order to identify the true ground state of the sys- 
tem it is necessary to take the atom-field interaction 
into account. Within the 4x4 degenerate subspace 
{\ip'ci),\ipa),\ipa),\tpi2)}^ the fuU Hamiltonian is given by 

/2lUc-2A -g -g \ 

-71-9 
-g 2uj,~2A 

\ -5 -71 2cu,-2aJ 

(20) 

This matrix can be diagonalized in closed form, yielding 
a unique ground state with energy Eg = 2ijJc — 2 A — VSg. 
The corresponding eigenstate is given by 

I'/'s) = 75 l^cl) + l^a) + i(IVa) + 1^:2))- (21) 

This approach constitutes a zeroth-order perturbation 
calculation in g for the case A = —A. The atom-cavity 
interaction has been included only to the extent that 
it lifts the ground state degeneracy. Nevertheless, the 
ground state (j2ip provides a good approximation when 
A ^ g. A comparison with the numerical solution of 
the ground state in the case A — lOg — —A is shown in 
Fig. [5] The two agree to within three percent. 

The primary lesson of this calculation is that the na- 
ture of the particles in the ground state of the coupled- 
cavity system is not necessarily determined by the atom- 
cavity interaction. The hopping term favors photonic 
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TABLE I: Characteristics of the four types of ground states in the coupled-cavity system. 



Phase 


ATVi 


Particles 


ANai 


Regime 


Insulator 





Atoms 
Polaritons 



> 


A/g<-l,^<lA| 
iA|/3<l,^<|A| 


Superfluid 


> 


Photons 
Polaritons 



> 


A<0,A>\A\>g;A>0,A>g 
A/g<^l,A^\A\ 



excitations over atomic. In the positive detuning regime, 
photons have lower energy than atomic excitations when 
A <^ g and so there is no competition between the in- 
teraction term and the hopping term. However, near 
A = the bare atom-cavity states are polaritonic. As 
A is increased the atomic component of the ground state 
is gradually eliminated, leaving a purely photonic super- 
fluid. At large negative detuning, the atom-cavity ground 
state is purely atomic. Near the superfluid boundary the 
hopping mixes the atomic and photonic components to 
create a polaritonlike superfluid. As the hopping is in- 
creased even further the ground state again reduces to a 
photonic superfluid. 

VI. CONCLUSIONS 

By examining the system over a wide range of values of 
the hopping and detuning parameters we have uncovered 
some unique features of the coupled-cavity model. The 
four types of states we have identified are summarized 
in Table HI They fall into two categories, analogous to 
the superfluid and insulator phases of the Bosc-Hubbard 
model, as indicated by the total excitation number vari- 



ance. Taking the atomic excitation number variance as 
an additional "order parameter" for the system allows the 
type of particles involved in the state to be determined. 
The insulator state may be either atomic or polaritonic, 
while the superfluid state may be photonic or polaritonic 
in nature. 

The different states arise from the bipartite nature of 
the system. Having both atoms and photons at each 
lattice site leads to competition between the atom-cavity 
interaction and the intercavity hopping. Both terms play 
a role in determining the phase of the system and the 
nature of the particles involved, although their relative 
importance depends on the particular parameter regime 
under consideration. Such richness of behavior suggests 
that the coupled-cavity model is more than just an analog 
for the Bose-Hubbard model and deserves further study 
in its own right. 
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